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(d  XI 


Notation 


=  inside  radius  of  tube. 

=  outer  radius  of  tube  walls, 
c  =  a  constant. 

d  =  minimum  distance  between  drops  in  spray  heat  exchanger, 
f  =  heat  transfer  coefficient  from  liquid  drops  to  salt  (1.42  watts/cm  -°C). 
g  =  acceleration  due  to  gravitational  force  of  the  earth  (980  cm/sec  ). 
h  =  power  transfer  coefficient  from  primary  coolant  to  tube  walls, 
k  =  thermal  conductivity  of  coolant  (Bi:  0.15  watts/cm-^C). 

I  =  length  of  each  tube. 

m  =  atoms  of  U  per  atom  of  Bi. 

n  =  number  of  drops  in  heat  exchanger, 

p  =  Laplace  transform  variable  replacing  t. 

Ap  =  pressure  drop. 

q  =  radial  distance  from  center  of  liquid  drop, 
r  =  radius  of  thermal  legs, 
s  =  iri  Vp/ ^ 

s-  =  roots  of  a  certain  transcendental  equation  (equation  (64)). 
t  =  time  taken  by  drop  to  fall  to  position  z  from  point  at  which  it  was  made, 
u  =  speed  of  coolant  in  heat  exchanger  tubes. 

V  =  speed  of  bismuth  in  drop-forming  nozzles. 

X  =  radius  of  header  at  point  y  along  its  axis, 
y  =  distance  along  header  from  beginning  of  flare -out. 

z  =  distance  below  plane  where  Bi  drops  are  formed,  at  which  the  drop  is  located, 
A  =  drop  radius, 

B  =  header  height. 

C  =  specific  heat  of  coolant  (Bi  coolant:  0.148  joules/g-  C). 

Cq  =  specific  heat  of  salt  (.92  joules/g-OC). 

C*  =  a  constant  =  0.44. 

D  “  CqMq/CM.  fk  •\ 

E  =  Young's  modulus  for  tube  walls  (Croloy  2  1/4:  29  x  10®  psi). 

G  =  number  of  holes  in  drop-forming  element  near  the  top  of  the  exchanger. 

H  =  power  flowing  through  each  exchanger  tube  wall. 

J  =  depth  of  Bi  pool. 

K  =  thermal  conductivity  of  tube  walls  (0.26  watts/®C-cm), 

L  =  total  length  of  thermal  legs,  i.e.  of  piping  connecting  reactor  with  exchanger. 
M  =  mass  rate  of  flow  of  the  primary  coolant,  viz,  Bi. 

Mq  =  mass  rate  of  flow  of  salt. 

N  =  number  of  tubes  in  heat  exchanger. 

P  =  power  generated  by  reactor  (450  Mw). 

Q  =  volumetric  flow  rate  of  bismuth. 

Qq  =  volumetric  flow  rate  of  salt, 

R  =  radius  of  the  heat  exchanger. 

S  =  stress  limit  under  operating  conditions  of  tube  walls  (Croloy  2  1/4:  8000  psi). 
T(A)  =  surface  temperature  of  drop. 

To(t)  =  salt  temperature  at  position  z  =  Ut. 

Xq(0)  =  temperature  of  salt  at  its  outlet. 

T(q,t)  =  drop  temperature  at  a  radial  distance  q  from  its  center  at  time  t. 

T(r,t)  =  average  temperature  of  drop. 
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A  SURVEY  OF  THREE  HEAT  EXCHANGERS  FOR  A  LIQUID 
METAL  FUEL  REACTOR 


Introduction 


Any  power -producing  reactor  must  be  cooled.  In  the  case  of  the  liquid  metal  fuel 
reactor  being  considered  by  Brookhaven  highly  radioactive  fuel,  U,  is  dissolved  in  a 
primary  coolant^ Bi,  of  low  cross  section.  The  primary  coolant  transfers  its  heat  by 
means  of  a  heat  exchanger  to  a  secondary  coolant,  which  operates  at  a  low  pressure. 
Finally,  the  secondary  coolant  transfers  its  heat  by  means  of  a  second  heat  exchanger 
to  a  third  thermodynamic  medium  operating  at  a  high  pressure.  At  the  present  stage 
in  the  art  of  Resigning  reactors  it  is  ill  advised  to  operate  the  secondary  coolant  at  a 
high  pressure  and  to  dispense  with  the  tertiary  coolant,  because  of  the  hazard  of  a  leak 
into  the  primary  coolant.  Also,  the  secondary  coolant  may  become  excessively  radio¬ 
active  to  be  used  as  the  final  thermodynamic  medium.  All  these  transfers  of  heat 
from  one  coolant  to  another  must  take  place  with  as  little  temperature  degradation  as 
possible  in  a  system  of  great  reliability( ^). 

The  primary  heat  exchanger  system  should  have  as  small  a  volume  as  is  conven¬ 
ient  in  other  respects  in  order  to  decrease  the  fuel  investment.  The  problem  of  how 
to  make  a  primary  exchanger  system  of  minimum  volume  or  at  least  of  small  volume 
is  the  concern  of  this  report.  We  do  not  consider  the  various  other  exchangers  which 
may  succeed  the  primary  one,  nor  do  we  consider  the  refinements  of  even  the  primary 
exchanger.  We  are  interested  rather  in  surveying  the  effects  of  the  many  quantities 
that  determine  the  exchanger  system  volume  and  general  configuration.  We  shall  try 
to  avoid  the  fine  points  which  would  have  to  be  considered  in  the  design  of  an  actual 
exchanger  system.  This  work  is  motivated  by  the  desire  to  get  a  better  estimate  of 
the  doubling  time  and  fuel  investment  in  our  liquid  metal  fuel  reactor  systems. 

In  this  report  we  shall  examine  three  heat  exchanger  systems  that  might  be  used 
in  conjunction  with  the  liquid  fuel  reactor  described  in  a  rudimentary  manner  in 
BNL-1u(2). 

1.  A  conventional  heat  exchanger  with  many  Croloy  Z  l/4  tubes  through  which 
the  primary  reactor  coolant,  the  bismuth -uranium  fuel  itself,  flows 
countercurrently  to  the  secondary  coolant.  The  primary  coolant  circulates 
by  reason  of  the  density  difference  in  that  coolant  between  the  hot  and  the 
cold  legs  connecting  the  heat  exchanger  and  the  reactor.  The  heat  ex¬ 
changer  is  located  high  enough  above  the  reactor  so  that  the  primary  cool¬ 
ant  circulates  at  the  requisite  rate.  We  call  this  manner  of  circulating 
the  coolant  "thermal  pumping." 


(^^The  author  hopes- that  some  day  all  this  temperature  degradation  and  complexity  can  be 
elimina-ted  by  some  one  magic  coolant.  Such  a  fluid  must  have  a  low  neutron  cross  sec¬ 
tion,  so  that  it  can  be  used  in  the  reactor  itself  and  will  not  become  significantly 
radioactive.  In  it  the  fuel  will  be  dissolved  or  suspended;  the  two  might  be  srpara- 
ted  easily  by  some  clever  method.  This  one,  prime  coolant  would  drive  turbines  di- 
rectlyf  the  fissionable  material  would  remain  behind  and  be  directly  returned  to  the 
reactor.  Heavy  water  or  he  1  i  um  unde  r  high  pressure  mi  ghtbe  a  candidate  fo  r  this  po  si  t  ion. 

Gurinsky,  and  others,  “Preliminary  Study  of  a  Uranium-Bismuth  Liquid  Fuel  Power 
Reactor, ”  BNL- 111,  (1951). 
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2.  The  second  system  is  the  same  as  that  in  scheme  1,  except  that  any 
thermal  pumping  is  augmented  by  a  mechanical  pump.  This  method  of 
circulating  the  coolant  is  called  mechanical  pumping. 

3.  The  third  scheme  uses  a  novel  heat  exchanger  in  which  the  very  dense 

bismuth  coolant  is  formed  into  drops  near  the  top  of  the  heat  exchanger. 
The  drops  are  allowed  to  fall  through  and  in  direct  contact  with  a  molten 
liquid  salt  of  low  density,  viz,,  LiCl-KCl.  This  salt  is  introduced  near 
the  bottom  of  the  heat  exchanger  cind  is  collected  at  the  top  above  the 
height  at  which  the  bismuth  drops  are  formed.  The  bismuth  drops  ag¬ 
glomerate  in  a  region  below  the  point  at  which  the  liquid  salt  is  introduced. 
The  salt,  in  addition  to  purifying  the  bismuth -uranium  fuel  of  all  fission 
products(3),  acts  as  the  secondary  coolant.  Thermal  pumping  appears  to 
be  out  of  the  question  with  this  scheme;  we  can  consider  only  mechanical 
pumping.  This  system  will  be  referred  to  as  the  "spray  exchanger 

system,"  The  other  type  of  heat  exchanger  we  shall  call  a  "tubular  heat 
exchanger."  We  shall  see  that  the  spray  type  of  heat  exchanger  has  many 
positive  virtues,  although  it  raises  several  problems. 

The  complete  power  system  may  be  divided  into  a  reactor,  a  heat  exchanger  which 
transfers  energy  from  the  primary,  radioactive  coolant  used  in  the  reactor  to  the 
secondary,  non-radioactive  coolant  used  as  a  thermodynamic  medium,  and  the  power - 
producing  unit  itself.  We  shall  consider  only  the  heat  exchanger,  the  importance  of 
various  factors  in  the  design  of  the  primary  heat  exchanger  system  and  their  effect 
on  the  investment  of  uranium.  We  may  characterize  the  reactor,  for  the  purposes  of 
this  report,  as  a  device  producing  a  certain  amount  of  power  at  a  specific  temperature 
differential  between  inlet  and  outlet,  and  requiring  a  certain  constant  pressure  differ¬ 
ential  to  drive  the  primary  coolant.  The  header  causes  a  loss  in  the  kinetic  energy  of 
the  coolant,  especially  in  the  flare -out. 

The  tubular  heat  exchanger  is  the  source  of  both  frictional  losses  and  kinetic 
losses  (especially  at  the  flare -out).  The  heat  energy  in  passing  through  the  walls  of 
the  heat  excheuiger  tube  is  degraded  in  temperature  both  as  a  result  of  the  gradient 
through  the  walls  and  of  the  drop  across  the  laminar  film  between  the  coolant  and  the 
tube  walls.  The  remaining  film  drop  between  the  outside  of  the  tubes  and  the  second¬ 
ary  coolant  is  left  to  the  consideration  of  the  engineer  designing  the  other  parts  of  the 
system  which  will  utilize  the  heat  energy.  Masses  of  uranium  are  invested  in  the  legs, 
the  header,  and  the  heat  exchanger.  The  investment  in  the  header  is  trivial,  however, 
compared  to  the  investment  elsewhere. 

For  reasons  of  convenience  we  choose  as  the  independent  variables  the  inner 
radius  of  the  heat  exchanger  tubes,  their  wall  thickness,  the  number  of  tubes  in  the 
heat  exchanger,  the  speed  in  the  thermal  legs,  the  temperature  drop  across  the  tube 
walls  and  the  laminar  film  between  the  walls  and  the  bismuth,  and  the  fractional 
kinetic  energy  loss  parameters.  From  these  arbitrarily  chosen  quantities  we  find  the 
radius  of  the  thermal  legs,  the  speed  in  the  heat  exchanger  tubes,  the  heat  transfer 
coefficient,  the  length  of  the  heat  exchanger  tubes,  the  volume  of  the  heat  exchanger 
system,  and  the  mass  of  uranium  invested  in  it.  In  the  case  of  thermal  pumping,  we 
calculate  the  length  of  the  thermal  legs.  In  the  case  of  mechanical  pumping,  this  leg 
length  is  an  independent  variable;  we  compute  the  pumping  power  and  pressure  differ¬ 
ential  the  pump  must  create. 


Bareis,  “A  Continuous  Fission  Product  Separation  Process,”  BNL-125,  (1951), 
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The  advantages  of  a  spray  exchanger  have  been  manifest  for  mamy  years:  the 
uranium  investment  is  relatively  small,  and  the  energy  is  available  at  but  a  small 
degradation  in  temperature.  In  this  exchanger  the  energy  is  degraded  in  temperature 
by  the  film  between  the  bismuth  drops  and  the  salt  and  by  the  gradient  in  the  drops 
themselves.  In  any  horizontal  plane  the  salt  is  considered  to  be  isothermal  because 
of  the  mixing  and  because  of  the  fact  that  the  drops  fall  uniformly  on  the  plane.  In 
this  exchanger,  the  bismuth  loses  all  its  kinetic  energy  in  being  formed  into  drops  at 
the  top  of  the  exchanger  and  also  suffers  a  modified  head  loss  in  falling  as  drops 
through  the  exchanger.  These  pressure  losses  preclude  the  use  of  thermal  pumping. 

Again  for  reasons  of  convenience  we  take  the  radius  of  the  drops,  the  speed  at 
which  the  salt  flows  through  the  heat  exchanger,  the  height  of  the  heat  exchanger  and 
the  ratio  of  the  mass  rates  of  flow  of  salt  to  bismuth  as  the  independent  variables. 
From  these  we  compute  the  speed  at  which  the  drops  fall,  the  difference  between  the 
inlet  and  outlet  salt  temperatures,  the  volumetric  and  mass  flow  rates  of  the  bismuth 
and  the  salt,  the  closest  spacing  of  the  grid  forming  the  drops,  the  number  of  drops, 
the  radius  of  the  heat  exchanger,  the  number  of  nozzles  needed  to  form  the  drops,  the 
investment  in  uranium  and  in  bismuth,  the  pressure  drop,  and  the  pumping  power. 

Several  of  the  various  parameters  seem  arbitrarily  chosen.  We  are  considering 
the  heat  exchanger  system  for  case  F  in  Table  A-1  of  BNL-111  and  make  the  as¬ 
sumptions  of  Appendix  1  of  that  report.  In  point  of  fact,  there  is  really  nothing  fixed 
about  the  power  developed  by  the  reactor,  or  the  pressure  drop  and  temperature 
differential  across  it.  They  could  be  very  different  from  the  numbers  selected  here. 

In  the  next  section  we  formulate  a  theory  for  tubular  heat  exchangers  and  for 
spray  heat  exchangers.  In  the  third  section  we  present  the  results  for  the  heat  ex¬ 
changers  considered  in  this  report. 


Formulation  of  the  Theory 


Tubular  Heat  Exchanger 

The  primary  coolant,  bismuth,  flows  as  the  heatant  inside  circular  tubes  over 
whose  external  surfaces  the  secondary  coolant  flows  countercurrently.  A  schematic 

diagram  is  shown  in  Figure  1.  If  the 
fluid  is  pumped  thermally,  then  the  heat 
uPAnFR— ^ HEADER  exchanger  must  be  high  enough  above 

- the  reactor  to  provide  a  pressure  differ 

^  ~  ential  between  the  hot  and  cold  legs  suf- 

^21-^  \  ^  LEG  ficient  to  circulate  the  bismuth  at  the 

LEG  \  / 

N  1  ^  rate  needed. 


HEAT  EXCHANGER 
—  THERMAL  LEGS  - 


REACTOR 


direction 
"of  coolant 
FLOW 


MECHANICAL 

PUMP 


Figure  1. 


The  mass  rate,  M,  of  flow  of  pri¬ 
mary  coolant,  bismuth,  is  fixed  by  the 
statement  that  the  heat  removed  from 
the  coolant  in  the  exchanger  is  equal  to 
the  heat  given  to  it  by  the  reactor: 

P  =  CMAT  ,  ( 


where  P  =  the  power  (450  Mw  in  our  case)  produced  by  the  reactor, 

C  =  the  specific  heat  of  the  coolant  (0,148  joules/g-°C), 

AT  =  the  temperature  difference  (ITO^C  in  our  case)  between  the  hot  and  cold 
thermal  legs,  each  of  which  is  assumed  isothermal,  i.e,,  no  heat  losses 
occur. 

The  temperature  drop,  ATj,  across  the  film  is  given  by  the  power,  H,  flowing 
through  each  tube  wall: 

H  =  ZiraJhATj  , 

=  p/n  . 


where  a  =  inside  radius  of  tube, 

I  =  length  of  each  tube, 

h  =  heat  transfer  coefficient  from  primary  coolant  to  tube  wall, 

N  =  number  of  tubes  in  the  heat  exchanger. 

The  temperature  drop,  AT^.  across  the  tube  wall  is  also  given  by  the  power  H  flowing 
through  each  tube: 

H  =  (4) 

in  b/a 

where  K  =  thermal  conductivity  of  tube  walls  (.26  watts/®C -cm), 
b  =  outer  radius  of  tube  walls, 

AT2  =  temperature  drop  across  tube  wall. 


10 


11 


The  heat  transfer  coefficient,  h,  is  found('^)  from  the  following  relation: 

2ha/k  =  7  +  0.025  (2Capu/k)^-®  ,  (5) 

where  k  =  thermal  conductivity  of  primary  coolant  (.15  watts/cm-^C), 
p  =  density  of  primary  coolant  (9.91  g/cm^)  , 
u  =  speed  of  coolant  in  heat  exchanger  tubes. 

The  speed,  u,  is  given  by  mass  conservation: 

Nira^up  =  irr^Vp  ,  (6) 


where  r  =  radius  of  thermal  legs, 

V  =  speed  of  coolant  in  thermal  legs. 

The  total  temperature  drop  ATj  +  AT^  which  can  be  tolerated  usually  limits  the 
amount  of  power  flowing  through  each  tube.  Under  some  conditions,  however,  the 
stress  the  tubes  can  endure  may  be  limiting.  For  thin  walls  the  stress,  S,  due  to  the 
thermal  gradient  is  given  by(^^ 


H  = 


4irS(l-t/)Kl 


t^E 


1  - 


2a2 


b^-a^ 


In  b/; 


(7) 


^  4TTS(l-v)Kl(b+a) 

77E(b-a) 

where  v  =  Poisson’s  ratio  (0.28), 

77  5=  thermal  coefficient  of  expansion  (16.  x  lO^^/^C), 

E  =  Young’s  modulus  (29  x  10^  Ibs/in^). 

We  need  to  take  into  account  the  fissionable  material  invested  in  the  two  headers. 
We  consider  a  header  shaped  so  that  the  momentum  of  the  primary  coolant  changes 
uniformly  with  distance.  (See  Figure  2.  We  shall  show  later  that  any  reasonable 
shape  of  the  header  gives  the  same  total  mass  in  the  heat  exchanger  systems.)  The 
shape  of  a  header  is  then  given  by 

y  =  c  [v  -  M/irpx^]  ,  (9) 

where  y  =  distcince  adong  axis  of  header  from  point  where  header  starts  flaring  out, 

X  =  radius  of  header  at  point  y, 
c  =  constant. 

The  value  c  may  be  determined  from  the  condition  that 


Lyon,  Liquid  Metals  Handbook.  Atomic  Energy  Commission,  Department  of  the  Navy, 
Washington,  1950,  Ch,  6,  p,  130  ff, 

(^)s,  Timoshenko,  Theory  of  Elasticity.  McGraw-Hill  Book  Co.,  1934,  para,  114. 
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dx 


=  tan  3.5°  =  ZcM/irpr^,  (10) 


V 


Figure  2. 


y=0 


which  is  the  condition  that  the  kinetic 
losses  be  minimized.  The  volume  of  the 
header  is  given  by  the  expression 


I  TTX^dy  , 
^0 


(11) 


where  B  =  header  height. 


Finally,  B  may  be  determined  from 
the  condition  that  the  header  cross- 
sectional  area  must  be  the  same  as  the 
area  of  all  the  tubes  inside  the  heat  ex¬ 
changer: 

trx^  =  Nira^  ,  (12) 

where 


xq  =  X  at  y  =  B, 


(We  neglect  the  small  header  volume  entailed  in  the  expansion  from  the  total  cross 
sectional  area  of  the  heat  exchanger  itself.)  The  mass  rate  of  flow,  M,  is  given  by 


M  =  irr^pV. 


(13) 


We  can  now  very  easily  find  the  volume  of  each  header  to  be 

0.19  r3  I  In  (v/ir a/r)  |  . 

The  volume  of  the  thermal  legs  is 


irr2L  , 


(14) 


(15) 


where  L  =  total  length  of  thermal  legs. 

The  volume  within  all  heat  exchanger  tubes  is: 

Nira^j . 


(16) 


The  total  volume  of  primary  coolant  in  the  heat  exchanger  system  is: 

2  X  0.19  r3  I  In  (<fW a/r)  I  +  irr^L  +  ira2Nj  . 


(17) 


The  mass  of  uranium  tied  up  in  the  heat  exchanger  system  is  easily  calculated 
to  be: 


(1 


_  m  [0.19  X  2r^  In  (i/Wsl/t)  I  +  irr^L  +  ira^Nl], 

,6  +  m)  I 


(18) 
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where  m  =  atomic  fraction  of  U  to  Bi* 

The  length,  L,  of  the  pipes  connecting  the  reactor  with  the  exchanger  must  be 
long  enough  to  go  from  one  to  the  other.  We  assume  a  shield  200  cm  thick  between 
the  reactor  and  heat  exchanger,  and  an  unblanketed  reactor  of  only  100  cm  radius  and 
of  150  cm  height.  We  assume  further  the  total  cross-sectional  area  of  the  exchanger 
to  be  four  times  as  great  as  the  area  within  the  exchanger  tubes  alone.  The  radius  of  the 
exchanger  is 


R  =  2  /n  a  , 

(19) 

The  minimum  active  pipe  length  is,  then,  from  Figure  1: 

L  >  2  [  2  \fW' a  +  100  +  20o] 

(20) 

^■4  a  +  600  . 

(21) 

The  minimum  passive  length  is 

1  1  -  150|  , 

(22) 

where  we  have  neglected  the  header  length.  Probably  the  pipes  should  be  rather  longer 
than  our  optimistic,  minimum  estimates. 


The  pressure  drop  due  to  the  frictional  losses  in  the  heat  exchanger  tubes  is 
given  by(6): 


Ap  =  1.52  X  10-6  jul.8/al.2  , 


(23) 


where  Ap  =  pressure  drop  in  psi,  I  and  a  are  expressed  in  cm.  u  is  in  cm/sec.  The 
pressure  drop  due  to  frictional  losses  in  the  thermal  legs  is  given  by 

Ap  =  1.52  X  10-6  i^yl.S/rl.2.  ,  (24) 


The  pressure  drop  due  to  kinetic  energy  losses  is  very  much  harder  to  compute. 
Such  losses  occur  primarily  in  the  flare-out  from  the  inlet  header  and  from  the  flare - 
out  of  each  heat  exchanger  tube  into  the  outlet  header.  We  circumvent  calculating  the 
kinetic  losses  in  characterizing  them  by  parameters  a  and  p;  p,  the  fraction  of  the 
kinetic  energy  per  unit  volume  of  the  thermal  legs  lost  in  the  headers  and  in  the  legs 
themselves,  and  a,  the  fraction  of  the  kinetic  energy  per  unit  volume  in  a  heat  ex¬ 
changer  tube  lost  in  the  entrance  and  exit  to  that  tube.  Were  the  kinetic  energy  of  the 
liquid  in  the  heat  exchanger  tubes  to  be  totally  lost  the  pressure  drop  would  be  pu^/2. 
The  total -pres sure  drop  due  to  kinetic  losses  is  then  characterized  by: 

Ap  =  7.25  X  10"^  |apu^  +  ppV^J  .  (25) 

The  pressure  generated  by  the  density  difference  of  the  coolant  between  the  hot  and 
cold  thermal  legs  is  given  approximately  by: 

Gilliland,  D,  Bareis,  G.  Feick,  in  the  Science  and  Engineering  of  Nuclear 
Power.  Addison  Wesley,  1949.  Vol.  II,  Ch.  10,  p.  173. 
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=  pL^ATg/2  ,  (2  6) 

=  ,00711  p^ATL  .  (27) 

where  ^  =  density  expansion  coefficient  of  the  primary  coolant,  Bi(l*22xl0  ^  g/cm^  -®C). 

This  expression  would  be  exact  if  L/2  were  the  altitude  difference  between  the 
horizontal  parts  of  the  pipes  joining  the  heat  exchanger  and  reactor  as  shown  in 
Figure  1,  If  these  horizontal  parts  are  short,  then  the  expression  will  be  a  good  ap¬ 
proximation.  These  horizontal  runs  should  be  as  short  as  is  reasonably  convenient 
to  keep  the  fuel  investment  low.  Let  us  assume  the  mechanical  pump  creates  a  pres¬ 
sure  difference  of  U  in  the  direction  of  the  pressure  created  by  the  thermal  legs.  After 
substitution  for  the  properties  of  bismuth,  the  pressure  balance  equation  then  becomes: 

n  +  1.47  X  10-3  L  =  1.52  x  lO'^  +  1.52  x  10'^  + 


7.19  X  10-3  (ou2  ^  pv^)  +  0.43  +  1.52  x  10-^  V^'®|  I  -  ISoj/r^'^  .  (28) 

We  are  now  in  a  position  to  set  up  the  procedure  for  study  of  the  heat  exchanger 
design.  The  choice  of  convenient,  independent  variables  is  governed  by  equation  (28). 

It  is  found  that  a,  b-a,  N,  V,  a,  p,  and  ATj  -H  AT^  are  the  most  convenient  set  of  in¬ 
dependent  variables.  In  the  case  of  mechanical  pumping  we  may  choose  L  in  addition 
provided  it  satisfies  the  inequality  (21).  Then  from  equations  (1)  and  (13)  r  may  be 
determined  as: 

r  =  .0357  (P/V)^/3  ,  (29) 


where  Pis  measured  in  watts  and  V  in  cm/sec.  We  combine  equation  (6)  and  (2  9)  to  get  u: 

u  =  1.27  X  10-3  (P/Na^)  .  (30) 


The  heat  transfer  coefficient  is  next  determined  from  equation  (5): 

541  u-^ 

h  =  --^-1.  .0203  — 8  , 
a  u‘° 


(31) 


where  a  is  taken  in  cm,  u  in  cm/ sec,  and  h  in  watts/cm^ -®C.  From  equations  (2), 
(3),  and  (4),  we  can  determine  the  length  of  heat  exchanger  tubes  needed  to  transmit 
the  power  desired: 


,  _  .159  (1  +  3.85  ha  in  b/a)  P 
(ATj  +  AT2)  Nah 

Nevertheless,  the  stress  limit  of  the  material  (Croloy  2.  l/4)  of  the  tube  walls  must 
not  be  exceeded.  From  equation  (7)  we  learn  that  the  I  found  from  equation  (33)  must 
satisfy  the  following  inequality: 

.0229  ^  (1  -  b/a)  .  (33) 

If  this  inequality  is  not  satisfied,  then  the  equality  of  (33)  is  used  to  determine  Z. 
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In  the  case  of  thermal  pumping,  the  total  length  of  the  thermal  legs  necessary  to 
circulate  the  coolant  at  the  requisite  rate  can  be  foundfrom  equation  (28)  by  putting  H  =  0: 


1.8 


1.8 


.43+  1,52x10"^  '  -150  +1.52x10"^  7.19x10"^  +  7.19x10-5  au^ 


L  = 


•a  L 

1.47x10'^  -  l.BZxlO-^-^^Y^ 


(34) 


However,  in  all  cases  the  independent  variables  must  be  adjusted  so  that  L  is  at  least 
as  great  as  the  minimum  necessary  length  given  by  equation  (21).  An  especially  con¬ 
venient  parameter  to  adjust  in  this  connection  is  V. 

In  the  case  of  mechanical  pumping,  the  length  of  the  pipes  connecting  the  reactor 
with  the  heat  exchanger  is  given  by  equation  (21).  Equation  (28)  reveals  that  the  pres¬ 
sure  which  must  be  created  by  the  mechanical  pump  is  given  by; 


n  =  0.43-1- 


1.52x10 


- -  1.47x10-3 


1.2 


1.8 


L+  1.52xlO-^-5^y^-l- 7.19x10-5  (cu^-l-pV' 


+  1.52x10 


.6  Y 


1.8 


1.2 


l~  150 


(35) 


The  pumping  power  is  the  work  done  per  unit  time,  W,  in  circulating  the  coolant  by 
the  mechanical  pump; 


W  =  Mn/p  ,  (36) 

=  2.76  X  10-5  OP  ,  (37) 

where  W  and  P  are  measured  in  watts  and  11  in  psi.  The  latter  equation  is  a  con¬ 
sequence  of  (1). 

In  either  case,  the  volume  of  the  heat  exchanger  system  is 

0.383  r3  I  in(>/ira/r)|  +  3,14  [Na^I  +  r2  (l +1  r  -  150  )],  (38) 

and  the  mass  of  uranium  is 

m  =  .00661  X  volume,  (39) 


Spray  Heat  Exchanger 

We  turn  now  to  the  quantitative  evaluation  of  the  spray  heat  exchanger.  Drops  of 
liquid  bismuth  fall  through  a  tank  full  of  KCl-LiCl  that  slowly  moves  upward.  See 
Figure  3  for  a  schematic  diagram. 

If  the  liquid  drops  of  Bi  are  assumed  to  fall  at  a  constant  speed  through  the  salt, 
the  calculation  is  greatly  simplified.  This  approximation  is  a  good  one.  The  speed  of 
fall,  Uf  is  given(^l  by  Newton*s  law; 


Perry,  Chemical  Engineer’s  Handbook,  Mc-Graw-Hill  Book  Co.,  1950.  3rd  ed,, 
p.  1017  ff. 
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Figure  3.  Spray  heat  exchanger. 


=  -Uq  +  176,0/A*,  (41) 

where  Uq  =  salt  velocity  (in  cm/ sec), 
pQ  =  salt  density  (1,60  g/cm^), 

C'  =  a  constant  (0.44), 
p  =  Bi  density  (9.91  g/cm^), 

A  =  Bi  drop  radius  (in  cm). 

The  heat  balance  equation  between  the 
Bi  drops  and  the  salt  states  that  the  heat 
transferred  through  the  laminar  film  sur¬ 
rounding  a  Bi  drop  per  unit  time  is  the 
same  as  the  heat  gained  by  the  salt  per 
unit  time: 


4nA2f  [T(A.t)-To(t)]  ^  ^  =  -CqPo  AzATo(t)  ,  (42) 


where  f  =  heat  transfer  coefficient  from  liquid  drops  to  salt  (1,42  watts/cm^-®C), 

T(A,t)  =  temperature  of  drop  at  its  surface  q  =  A  at  time  t, 

X  =  active  length  of  spray  heat  exchanger, 
n  =  number  of  drops  in  heat  exchanger, 

Cq  =  specific  heat  of  salt  (0,92  joules/g-®C), 

Tq  =  salt  temperature, 

z  =  distance  below  horizontal  plane  where  drops  are  formed, 


R  =  radius  of  the  heat  exchanger  . 


Let  the  initial  uniform  drop  temperature  be  taken  as  zero.  The  drop  temperature 
at  any  later  position  is  given^®)  by: 

CO 

X(q,t)  =  ^  (-)j  j  1  d^  exp  {kj^Tr^(  ^ -t)/pCA^)T(A,  O  sin  (irjq/A)  ,  (43) 

pCAq 

where  T(q,t)  =  drop  temperature  at  point  q  from  its  center  at  a  time  t, 

t  =  time  taken  by  drop  to  fall  to  a  position  z  from  point  at  which  it  was 
made , 

k  =  thermal  conductivity  of  Bi  (0,15  watts/ cm-®C), 


Note  that 

z  =  Ut  ,  (44) 

To  avoid  making  mistakes  it  is  most  essential  that  equation  (43)  be  converted  to  a  uni¬ 
formly  convergent  series.  Upon  integrating  by  parts,  we  get  such  a  series: 

T(q.t)-T(A.t)  =-^  ^  ^/o  ^  -tVpCA^)  sin  (irjq/A)  ,  (45) 

j=l 


^®)h.S.  Carslaw  and  JT.C.  Jaeger.  Conduction  of  Heat  in  Solids*  Oxford.  1947.  para.  91* 
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where  we  have  used  the  familiar  relation: 

2  sin  (jirq/A)  =  -Trq/2A  • 

j=l  J 

The  heat  in  the  drops  in  excess  of  that  in  them  at  zero  of  temperature  is 

4TrpC  I  T(q,t)q^dq 


=  ^-trpCA^T(A.t) 

3  ‘fT 


.CA^  ”  _i.r* 
—  j=i 


exp  (kj2ir2(  I  -t)/pCA2)  .  (48) 


The  heat  energy  leaving  a  drop  per  unit  time  must  equal  the  heat  absorbed  by  the  salt 
per  unit  time: 

-4xpC  A  p  T(q.t)q2dq  =  .  (49] 

Jq 


where  Mq  =  mass  rate  of  salt  flow, 


4  on 


Mq  =  PqUoIt^J^  -  — ttA 


3  \  /  . 


This  last  relationship  may  be  immediately  integrated.  In  conjunction  with  (48)  we  get: 


[W  -  W]  =  -4iTpC  T(q.t)q2dq 

Jo 


=  --|-TrA^pCT(A,t)  +  2  ^  r  d^  exp  (kj2Tr2(  4 -t)/pCA2) - (  — (52) 

3  "  j=l  J  Jo 


where  Tq{0)  =  salt  temperature  at  exit. 

Equation  (42)  may  be  solved  for  T(A,t): 

T(A  t)  -  T  (t)  -  Ao  .  (< 

T(A.t)-To(t)  ^^2nfu  dt  ‘ 

By  substitution  of  (53)  into  (52)  we  can  get  a  simple  integral  equation  in  TQ(t)  alone; 


J 


-/cn^t 


”  [  A  J_+  ^ApCpMo  n  p  ^^gKn2(C-t) 
j=lU2K  irZfM  j^J  Jo  dC^ 
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where 


D  =  CoMo/CM, 

M  =  mass  rate  of  Bi  flow, 
4  o  n 

=  ^TrA^p-j^U, 

K  -  kirV  pCA^  . 


(55) 

(56) 

(57) 


The  integral  equation,  (54),  may  be  solved  by  the  Laplace  transform  method;  (p 
is  the  variable  conjugate  to  t*) 

[pT5<p)-'fo(»)]  *  ik  * 


CO  r  dTo{0)1  r  1  ]  [  6  ,  ZApCpMo] 

.2  [p2To{p)-pTo(0)  -  J  [^z\  [j4„2^  +  ,2fMj2  J  '  °  ' 

3 


(58) 


We  use 


TT  K 


j=l  j^(P+^j^)  ^P^  2p^  p 


coth 


(59) 


TT^K  _  K‘ 


j_l  j4(p+/<j^)  90p  6p2  2p 


(60) 


to  evaluate  several  series.  The  Laplace  transform  of  the  salt  temperature  may  now 
be  found; 


To(p)  = 


D  1 

[l-k/fA] 

1  DTo(0)/J'cothv^ 

fz 

K 

D[l-k/fA] 

3; 

^  --J 

'P  TT  V 

f — coth  itJ 
P  ^ 

rjcothj 

K 

(61) 


! 


upon  using  the  fact  that  T(a,0)  =  0  in  equation  (63)  to  evaluate  dTQ(0)/dt  in  terms  of 
Tg(0).  The  salt  temperature  is  found  from  the  inverse  of  (61).  Let 

s  =  iri  ^p/ K  .  (^^) 


If  D  /  1; 


Tn(t)  3  _2. 
To(0)  D-1 


00 

6D  s 

j 


KS^t/-n* 

J  ' 


i=l  [(3.2D  +  ^  s2)  (3  +  g  s|)  -  3D  +  d2  (l  +  sj 


(63) 


where  s.  are  the  nonzero  solutions  of  • 


i 


1 
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(64) 


Sj  cot  Sj  =  1  - 


fA  ^  Ds^ 


If  D  =  1,  there  exists  a  double  pole  in  (61 )  at  p  =  0,  The  salt  temperature  is  then 


given  by: 


-  ^s 


Tn(t) 

T 


ftA2 


The  surface  temperature  of  the  metal  is  found  from  (63)  and  (53).  If  D  ^  1: 


T(A,t)  ^  _D_ 
To(0)  D-1 


6D 


^  ^  3fA 


^  / 
ijt/ir' 


j=l 


^3-20  +  —  s?\  ^3+—  s2\  -  3D  +  D2  (l+k/fA)s? 
\  iA  })\  iA  1/  J 


(65) 


(66) 


Equation  (45)  may  be  used  to  find  T(q,t)if  desired.  However,  the  average  temperature, 
T(q,t),  is  of  more  interest.  The  use  of  equation  (51)  leads  to  considerably  simpler  re¬ 
sults  than  the  more  straightforward  method,  viz,: 


y  wA^pC 


(67) 


=  D(To(t)  -  To(0))  ,  (68) 

into  which  (63)  may  be  substituted  if  desired. 

Let  us  turn  for  a  moment  to  a  consideration  of  the  roots  of  (64).  First,  let  us 
look  at  a  rough  plot  of  the  right  hand  side  and  left  hand  side  of  (64)  in  Figures  4a  and 
4b,  The  following  points  are  to  be  noted: 

1,  The  s  cot  s  lines  go  through  the  odd  half  interger  points  with  increasing 
steepness  as  s  gets  larger, 

2,  Both  the  right  side  and  left  side  of  the  equation  start  out  at  s  =  0  with  unit 
value  and  zero  slope. 

3,  The  1-1/  function  approaches  1  -  asymptotically. 

4,  The  first  purely  real  root  occurs  at  s  >  it  if  D  <  1  but  ats<iTif  D  >1, 

5,  The  roots  are  spaced  at  least  it  apart,  but  this  spacing  approaches  exactly 
IT  as  s  increases  and  the  roots  become  very  nearly  a  half  odd  interger  xn 
as  s  gets  larger. 

Let  us  look  for  pure  imaginary  roots  now  in  Figures  4c  and  4d.  If  D  >  1  it  is  easy 
to  see  from  the  second  derivative  of  both  sides  of  equation  (64)  that  there  are  no  purely 
imaginary  roots.  If  ID  <  1,  by  the  same  method,  we  can  see  there  is  exactly  one  purely 
imaginary  root. 
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As  yet  we  have  not  considered  whether  there  may  be  complex  roots.  If  equation 
(64)  be  expanded  in  a  power  series,  it  can  be  reduced  to 


We  now  note  that  the  zeros  of 

<™' 

are  all  real  for  the  discriminant  is  always  positive.  If  D  <  1,  one  root  is  positive  and 
the  other  is  negative.  If  D  >  1,  both  roots  are  negative.  From  a  slight  generalization 
of  several  theorems  derived  by  Laguerre  relating  to  the  zeros  of  integral  functions, 
we  can  conclude(9)  that  in  the  p  plane  the  roots  of  the  transcendental  equation  are  all 
real.  Further,  in  the  p  plane  if  D  <  1,  one  root  is  positive  and  all  the  others  are 
negative,  and  if  D  >  1,  all  roots  are  negative.  The  summation  in  equation  (63)  is  to  be 
carried  over  all  roots  of  (64)  in  the  p  plane.  This  statement  means  that  we  sum  (63) 
over  all  non-zero,  negative  roots  Sj  of  (64).  Because  (64)  is  even  between  positive 
and  negative  values,  we  may  alternatively  sum  (63)  over  just  the  positive  solutions  of 
(64).  If  D  <1,  we  have,  of  course,  to  sum  (63)  over  the  one  imaginary  root  of  (64), 
also. 

We  may  now  evaluate  all  other  quantities  of  interest  in  sequence:  U  is  found  from 
equation  (36).  The  mass  flow  rate,  M,  of  bismuth  is: 

M  =  P/C6T  ,  (71) 

where  5T  =  the  differential  between  the  inlet  and  outlet  bismuth  temperatures  (ITO^C), 

6T  =  T(q,t)  .  (72) 

The  volumetric  flow  rate,  Q,  of  the  bismuth  is: 

Q  =  M/p  .  (73) 

The  mass  flow  rate,  Mq,  of  the  salt  is 

Mo  =  D(C/Co)M  ,  (74) 

The  volumetric  flow  rate  of  the  salt  is 

Qo  =  “o/po  •  (75) 

The  minimum  distance,  d,  between  drops  is 


The  number  of  drops,  n,  is 

C^^E.C.  Titchmarsh,  The  Theory  of  Functions,  2nd  edition,  para*  8*61  -  8,63, 
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n  =  3X.Q/ 4TrA^U  • 


(77) 


The  radius,  R,  of  the  heat  exchanger,  based  on  the  assumption  that  the  drops  are 
spaced  on  a  square  lattice  the  minimum  distance  apart,  is: 

R  =  (nd3Ax)l/2  .  (78) 

The  number  of  holes,  G,  in  the  top  plate  is: 

G  =  itR2/d2  .  (79) 

Again  if  we  assume  a  shield  of  200  cm  thickness,  a  reactor  of  100  cm  radius  and 
of  200  cm  height,  the  total  length  of  pipe  connecting  the  best  exchanger  and  reactor 
cannot  be  less  than 

L  =  600  +  2R  +  I  J  +  \  -  200  |  ,  (80) 

where  J  =  depth  of  the  bismuth  pool.  The  mass,  (i.  of  bismuth  in  the  spray  heat  ex- 
changer  is 

^  =  Mt  ,  (81) 

and  the  mass,  X,  of  uranium  in  the  exhanger  is 

X  =  (233/209)iim  .  (82) 

The  mass  of  bismuth  in  the  pool  where  the  drops  are  agglomerated  at  the  bottom 
of  the  heat  exchanger  is  ■jrJR2. 

Again,  the  mass  of  bismuth  in  the  pipes  connecting  the  reactor  and  heat  exchanger 
is  ■irr2Li. 

The  masses  of  uranium  in  the  pool  and  in  the  pipes  are  respectively  [l8,3/(1.60+ 
m)]  miTJR2,  and  [l8.3/(1.60+m)]  mirr2L. 

If  we  assume  the  drops  are  of  the  same  diameter  as  the  nozzles  or  hole  which 
formed  them,  the  speed,  v,  in  the  nozzle  is  approximately 

V  =  Q/GirA2  .  (83) 

Consequently,  the  pressure,  0,  that  must  be  created  by  a  mechanical  pump  if  we 
neglect  any  thermal  pumping  is: 

Ft  1 

n  =  1.45  X  10-8  j^Ypv2  +  gX(p-pQ)  +  ppv2/2  |  +  1.52  x  10  ^1.2  ^  »  (84) 

where  H  is  in  psi,  but  all  other  quantities  are  in  the  cgs  system  of  units.  Finally, 


w  =  no . 


(85) 


Results 


The  results  of  our  numerical  work  are  presented  in  graphical  form  in  Figures  5 
through  34.  The  results  may  be  broken  down  into: 

1,  any  of  the  heat  exchanger  systems, 

E.  any  tubular  exchanger  system, 

3.  only  tubular  exchangers  with  thermal  pumping, 

4.  only  tubular  exchangers  with  mechanical  pumping, 

5.  only  spray  exchangers. 

We  summarize  our  results  below. 


Any  Exchanger  System 

Radius  of  Piping  Connecting  Reactor  with  Exchanger:  Figure  6. 

1,  Decreases  as  (speed  in  piping)" It  is  desirable  to  use  as  high  a  speed  as 
is  convenient  to  cut  down  the  fuel  investment. 


Any  Tubular  Exchanger 

Speed  in  Heat  Exchanger  Tubes:  Figures  6a,  b. 

1,  Decreases  rapidly  as  (number  of  tubes) 

2.  Decreases  rapidly  as  (radius  of  exchanger  tubes) 

Power  Transfer  Coefficient:  Figures  7a,  b,  c. 

1,  Decreases  rapidly  as  tube  radius  increases.  Turbulence  decreases. 

2,  Increases  as  speed  in  tubes  increases,  but  more  slowly  than  linearly.  ^ 

3,  Decreases  as  number  of  tubes  increases,  but  slower  than  (number  of  tubes)"  . 
Speed  in  tubes  decreases. 

Length  of  Heat  Exchanger  Tubes:  Figures  8a,  b. 

1.  Decreases  more  slowly  than  (number  of  tubes) if  the  temperature  drop 
across  inner  laminar  film  +  wall  is  limiting;  decreases  as  (number  of  tubes) 
if  the  stress  the  tubes  can  endure  limits  (large  temperature  drops). 

2.  Decreases  somewhat  as  radius  of  tubes  increases  if  the  stress  limits.  Only 
I  slight  decrease  if  temperature  drop  limits. 

3.  Decreases  as  walls  become  thinner. 

Tube  Wall/inner  Laminar  Film  Temperature  Drop:  Figures  9a,  b. 

1.  Decreases  moderately  as  tubes  become  larger. 

2.  Decreases  as  number  of  tubes  increases. 

3.  Decreases  as  walls  become  thinner. 


Tubular  Exchanger  Thermal  Pumping 

Total  Active  Length  of  Pipes  Connecting  Reactor  with  Exchanger:  Figures  10a,  b,  c,d. 
1.  Decreases  extremely  rapidly  as  number  of  tubes  increases  up  to  10,000  and 
then  levels  off. 
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2.  Decreases  very  rapidly  as  radius  of  tubes  increases  up  to  1.0  cm,  then  levels 
off. 

3.  Increases  very  rapidly  as  speed  in  piping  increases,  approaching  an  asymptote 
at  480  cm/sec  =16.  ft/sec  beyond  which  friction  exceeds  density  difference. 

4.  Thinner  walls,  higher  temperature  drops  and  less  kinetic  energy  loss  are  of 
secondary  importance. 

Volume  of.  Bismuth  Weight  in.  Uranium  Investment  in  Whole  Heat  Exchanger  System: 

Figures  11a,  b,  c,  d. 

1.  Decrease  extremely  rapidly  as  number  of  tubes  increases  up  to  10,000  be- 
ca^use  piping  length  decreases  rapidly  and  then  levels  off, 

2.  Decrease  very  rapidly  with  tube  radius  to  a  minimum  at  1,0  cm  radius,  then 
increases  substantially.  At  first,  length  of  piping  decreases  very  rapidly,  but 
then  surface -to-volume  ratio  decreases, 

3.  From  a  broad  minimum  at  300  cm/sec  increase  rapidly  to  oo  at  480  cm/sec 
as  speed  in  piping  increases.  Low  speeds  imply  a  large  pipe  diameter  and 
high  speeds  high  friction;  in  between  there  is  an  optimum. 

4.  Increase  as  tubes  grow  fatter,  temperature  drop  grows  smaller,  and  as  kinetic 
energy  losses  increase, 

5.  Header  volume  is  negligible. 


Volume  of.  Bismuth  Weight  in,  and  Uranium  Investment  in  Heat  Exchanger  Proper: 

Figures  12a,  b, 

1,  Increases  rapidly  as  radius  of  tubes  increases  because  of  lower  surface-to 
volume  ratio, 

2,  Increase  moderately  as  number  of  tubes  increases  if  temperature  limits,  but 
remain  constant  if  stress  limits, 

3.  Decrease  as  walls  get  thinner  and  as  temperature  drop  grows  larger, 

4.  This  volume  is  important. 

Volume  of.  Bismuth  Weight  in,  and  Uranium  Investment  in  Piping  Connecting  Reactor 

with  Exchanger:  Figures  13a,  b,  c,  d, 

1.  Decrease  extremely  rapidly  as  number  of  tubes  in  exchanger  increases, 
leveling  off  at  10,000,  because  piping  length  decreases. 

2.  Decrease  very  rapidly  as  tube  radius  increases,  because  piping  length  de¬ 
creases,  and  level  out  at  1,0  cm  radius. 

3.  From  a  minimum  near  300  cm/sec,  increase  rapidly  to  oo  at  480  cm/sec  be- 
cause  of  friction, 

4.  Decrease  as  temperature  drop  increases,  as  walls  become  thinner,  and  as 
kinetic  losses  are  lowered, 

5.  This  volume  is  important. 


Tubular  Exchanger  Mechanical  Pumping 

Length  of  Piping  Connecting  Reactor  with  Exchanger:  Figures  14a.  b. 

1.  Decreases  rapidly  with  increase  in  number  of  tubes  to  a  broad  minimum  at 
10,000  tubes. 

2.  Longer  the  fatter  the  tube  walls,  expepially  for  small  radii, 

3.  Shorter  the  larger  the  temperature  drop,  especially  for  small  tubes. 

4.  Shallow  minimum  as  tube  diameter  is  altered. 

5.  This  is  the  weakest  part  of  the  calculation  because  one  may  want  to  add  many 
things  between  the  reactor,  its  shield  and  the  exchanger,  and  because  of  the 
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crudity  of  our  estimate  of  exchanger  radius. 

Total  Pressure  Drop:  Figures  15a,  b,  c,  d,  e. 

lo  Decreases  extremely  rapidly  as  number  of  tubes  increases,  levelling  out  at 
10,000  tubes. 

2,  Decreases  very  rapidly  as  tubes  grow  in  size  to  a  broad  minimum  at  1,5  Cm 
radius. 

3,  Increases  rapidly  with  speed  in  the  piping, 

4.  Decreases  a  little  for  thinner  walls,  especially  for  tubes  of  small  radius,  for 
higher  temperature  drops,  for  shorter  piping,  and  for  lower  kinetic  energy 
losses. 

5.  It  would  appear  that  pressure  limitations  are  reached  before  pumping  power 
limitations. 

Pressure  Drop  Caused  by  Exchanger  Tube  Friction:  Figures  l6a,  b,  c. 

1,  Decreases  extremely  rapidly  as  number  of  tubes  increases. 

2,  Decreases  very  rapidly  as  tubes  get  larger, 

3,  Decreases  a  little  for  thinner  tubes  and  higher  temperature  drops. 

4,  This  pressure  drop  is  important. 

Pressure  Drop  Caused  by  Kinetic  Energy  Loss  in  Heat  Exchanger  Tubes:  Figures 
17a,  b. 

1,  Decreases  extremely  rapidly  as  number  of  tubes  increases. 

2,  Decreases  very  rapidly  as  radius  of  tubes  increases, 

3.  Pressure  drop  is  of  minor  importance. 

Pressure  Drop  Caused  by  Friction  in  Active  Part  of  Piping  Connecting  Reactor  with 
Exchanger:  Figures  18a,  b,  c. 

1,  Increases  rapidly  as  speed  in  piping  increases, 

2.  Increases  as  heat  exchanger  tube  radius  increases  because  of  larger  diameter 
exchanger.  Increases  for  pressure  generated. 

3*  Increases  a  little  as  the  number  of  exchanger  tubes  increases.  The  pressure 
generated  behaves  similarly, 

4.  This  pressure  drop  is  of  minor  importance. 

Pressure  Drop  Caused  by  Kinetic  Energy  Loss  in  Piping  Connecting  Reactor  with 
Exchanger:  Figure  19. 

1.  Increases  rapidly  as  (speed  in  piping)^, 

2,  This  pressure  drop  is  a  major  contribution  to  the  total  drop. 

Pressure  Drop  Caused  by  Friction  in  Passive  Part  of  Piping  Connecting  Reactor  with 
Exchanger:  Figures  20a,  b,  c,  d, 

1.  Decreases  rapidly  as  number  of  tubes  increases, 

2.  Increases  rapidly  as  speed  in  piping  increases, 

3.  Decreases  moderately  as  size  of  exchanger  tubes  increases, 

4.  Increases  moderately  as  length  of  piping  increases. 

5.  Decreases  as  tubes  becomes  thinner  and  as  temperature  drop  increases. 

6.  This  pressure  drop  is  minor. 

Volume  of,  Bismuth  Weight  in,  and  Uranium  Investment  in  Heat  Exchanger  System; 
Figures  21a,  b,  c,  d. 

1.  Decrease  rapidly  as  number  of  tubes  in  increased  til  a  minimum  is  reached 
at  about  10,000  tubes. 
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2,  Increase  rapidly  as  the  radius  of  the  exchanger  tubes  increases  because  of 
the  lower  surface -to -volume  ratio.  The  lower  limit  on  tube  size  is  fixed  by 
limits  on  the  frictional  pressure  drop  and  by  requirements  of  rigidity. 

3,  Increase  as  tubes  walls  become  fatter  and  as  temperature  drop  decreases. 

4,  Increase  linearly  with  length  of  piping. 

5,  Decrease  a  little  as  speed  in  piping  increases. 

Volume  of.  Bismuth  Weight  in,  and  Uranium  Investment  in  Piping  Connecting  Reactor 

with  Exchanger:  '  Figures  22a,  b,  c. 

1.  Decrease  as  (speed  in  piping)* 

2.  Decrease  somewhat  as  number  of  exchanger  tubes  increases  and  reach  a  broad 
minimum  at  about  10,000  tubes. 

3.  Go  through  a  shallow  minimum  at  roughly  0.75  cm  as  the  tube  radius  is  in¬ 
creased. 

4.  Increase  as  tube  walls  become  thicker  and  as  temperature  drop  decreases. 

Volume  of  Bismuth  Weight  in,  and  Uranium  Investment  in  Heat  Exchanger  Proper: 

Figures  23a,  b, 

1,  Increase  rapidly  as  the  tube  radius  increases. 

2,  Increase  somewhat  as  the  number  of  exchanger  tubes  increases  if  the  tempera¬ 
ture  drop  limits,  but  remains  constant  if  the  stress  limits. 

3,  Increase  somewhat  as  the  wall  thickness  increases  and  as  the  temperature 
drop  increases. 

Volume  of.  Bismuth  Weight  in,  and  Uranium  Investment  in  Headers:  Figures  24a,  b,  c. 

1.  This  volume  is  strictly  of  minor  importance. 


Spray  Exchanger  Mechanical  Pumping 


Before  summarizing  the  results,  it  should  be  stated  that  a  weakness  in  our 
calculations  lies  in  the  fact  that  the  power  transfer  coefficient  has  conservatively 
been  assumed  to  be  1.42  watts/ cm^-OC,  instead  of  having  been  calculated  from  the 
relative  drop  and  salt  speed,  the  drop  size,  and  the  properties  of  bismuth  and  LiCl- 
KCl.  Measurements  will  be  started  to  rectify  this  situation. 

Temperature  Drop  From  Bismuth  to  Salt,  Hot  Pipes:  Figures  25a,  b,  c,  d. 

1.  Decreases  as  the  time  of  fall  of  bismuth  drops  increases. 

2.  Decreases  rapidly  as  the  heat  absorption  ratio  decreases.  The  heat  absorption 
ratio  is  the  ratio  of  heat  absorbed  by  salt  per  degree  to  the  heat  absorbed  by 
bismuth  per  degree. 

3.  Decreases  rapidly  as  the  drop  size  decreases  especially  for  low  heat  ab¬ 
sorption  ratios. 

Inlet  Salt  Temperature:  Figures  26a,  b,  c,  d,  e. 

1,  Increases  as  the  time  of  drop  fall  increases. 

2,  Increases  as  the  drop  size  decreases  for  small  times  but  not  for  large  time. 

3,  Increases  as  heat  absorption  ratio  increases. 

4,  Large  drops  come  into  thermal  equilibrium  after  2  sec  but  small  drops  come 
into  thermal  equilibrium  essentially  immediately. 
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Height  of  Active  Part  of  Heat  Exchanger:  Figures  27a,  b,  c. 

1.  Increases  moderately  as  time  of  fall  increases. 

2.  Increases  as  the  drop  size  increases. 

3.  Decreases  as  salt  speed  increases. 

4.  Note  the  active  part  of  the  exchanger  is  very  shallow. 

Minimum  Radius  of  Heat  Exchanger;  Figures  28a,  b.  .  .  , 

1.  Decreases  a  little  as  drop  size  increases,  as  heat  absorption  ratio  decreases, 
and,  for  the  larger  drops,  as  salt  speed  increases. 

Number  of  Holes  in  Distribution  Plate;  Figures  29a,  b. 

1.  Decreases  rapidly  as  the  drop  size  increases.  The  maximum  number  of  holes 

fixes  the  smallest  drop  size  allowable. 

2,  Is  almost  constant  as  the  salt  speed  changes  and  as  the  heat  absorption  ratio  changes. 

Total  Pressure  Drop;  Figures  30a,  b,  c,  d,  e,  f. 

1.  Increases  rapidly  as  the  speed  in  the  piping  increases. 

2.  Increases  somewhat  as  the  drop  size  increases. 

3.  Increases  as  time  of  fall  increases. 

4.  Decreases  somewhat  as  the  salt  speed  increases. 

5.  Is  almost  constant  as  the  heat  absorption  ratio  changes. 

Volume  of.  Bismuth  Weight  in,  and  Uranium  Investment  in  Heat  Exchanger;  Figures 
31a,  b,  c,  d,  e,  f,  g,  h. 

1.  Increase  somewhat  as  drop  sizes  decrease. 

2.  Increase  somewhat  as  time  of  fall  of  drops  increases. 

3.  Increase  a  little  as  heat  absorption  ratio  increases. 

4.  Is  almost  independent  of  salt  speed, 

5.  Decrease  a  little  as  speed  in  piping  increases. 

Volume  of,  Bismuth  Weight  in,  and  Uranium  Investment  in  Piping  Connecting  Reactor 
with  Exchanger;  Figures  32a,  b,  c,  d, 

1.  Decrease  a  little  as  drop  size  increases. 

2.  Increase  but  little  as  the  heat  absorption  ratio  increases, 

3.  Is  almost  independent  of  time  of  drop  fall. 

Volume  of,  Bismuth  Weight  in,  and  Uranium  Investment  in  Pool;  Figures  33a,  b. 

1,  Decrease  somewhat  as  drop  radius  increases. 

2,  Decrease  as  salt  speed  increases. 

3,  Decrease  a  little  as  heat  absorption  ratio  increases. 

Volume  of.  Bismuth  Weight  in,  and  Uranium  Investment  in  Active  Part  of  Heat  Ex- 
changer;  Figure  34. 

1,  Increase  markedly  with  time  of  drop  fall. 


Summary  of  Results 


We  summarize  our  results  in  Table  1  by  comparing  a  reasonable  design  of  each 
of  the  three  exchanger  systems.  We  note  the  superiority  of  the  spray  exchanger  over 
the  tubular  exchanger  with  mechanical  pumping  and  of  the  latter  over  the  tubular 
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Passive  part  .10  psi  ,23  psi 

Kinetic  (20%  energy  loss)  1,3  psi  2.30  psi  Kinetic  (20%  energy  loss)  3.6  psi 
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exchanger  with  thermal  pumping.  Because  the  length  of  piping  is  very  uncertain  we 
quote  the  net  pressure  created  per  unit  length  and  the  investment  per  unit  length. 

These  quantities  allow  us  to  take  into  account  bends,  which  were  omitted  from  con¬ 
sideration,  the  effects  of  inserting  various  components  which  were  not  considered, 
and  consequences  of  changing  the  reactor  pressure  drop  from  the  very  hypothetical 
,43  psi  to  some  more  realistic  value. 

It  may  be  noted  in  the  spray  exchanger  that  if  the  temperature  drop  from  the 
bismuth  to  salt  is  increased  to  21°C  the  time  of  fall  is  cut  to  1  sec  and  the  heat  ex¬ 
changer  investment  is  cut  in  half.  Remaking  the  pool  from  a  right  cylinder  into  a 
trumpet  will  materially  decrease  the  pool  investment.  The  spray  exchanger  then 
looks  even  more  advantageous,  (It  must  be  remembered  that  one  film  +  wall  tempera¬ 
ture  drop  with  the  tubular  exchangers  is  35°C,  whereas  the  total  drop  may  be  about 
50^0.)  There  are  few  other  parameters  which  may  as  effectively  lower  the  investment 
(although  a  higher  value  for  the  power  transfer  coefficient  of  the  drops  will  also  de¬ 
crease  the  spray  exchanger  investment).  The  disadvantages  of  the  spray  exchanger 
system  are:  (1)  difficulties  with  emergency  cooling,  (2)  difficult  corrosion  problems, 
(3)  possible  admixing  of  bismuth  and  salt,  (4)  somewhat  higher  pumping  powers  (not 
serious),  (5)  number  of  holes  in  the  distribution  plate. 

Kinetic  energy  losses  in  the  heat  exchauigers  proper  are  of  negligible  importance. 
Header  volumes  are  of  negligible  importance. 
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Any  Tubular  Exchanger 


vs  radius  of  heat  exchanger  tubes. 


WALL  THICKNESS, 

DROP  ACROSS  INNER  FILM  AND  TUBE  WALL 


Tubular  Exchanger 
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INNER  DIAMETER  OF  TUBES  (INI  Figure  9b.  Tube  wall  film  temperature 

0.25  I  0.5  .75  f.O  I  1.25  1.5 1  175  2.0  drop  VS  number  of  heat  exchanger  tubes. 


Tubular  Exchanger:  Thermal  Pumping 
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10a.  Total  active  length  of  pipes  connecting  Figure  10b.  Total  active  length  of  pipes  connecting 

with  heat  exchanger  vs  number  of  tubes.  reactor  with  heat  exchanger  vs  inner  radius  of  tubes. 


Tubular  Exchanger:  Thermal  Pumping 
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Figure  10c.  Total  active  length  of  pipes  connecting  Figure  lOd.  Total  active  length  of  pipes  connecting 

reactor  with  heat  exchanger  vs  inner  radius  of  tubes,  reactor  with  heat  exchanger  vs  speed  in  pipes. 
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Figure  12a.  Mass  of  uranium  in  and  volume  of 
heat  exchanger  proper  vs  inner  radius  of  heat 
exchanger  tubes.  BNL  Log  No.  D-2072, 


Tubular  Exchanger:  Thermal  Pumping 


Tubular  Exchanger:  Thermal  Pumping 
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Tubular  Exchanger:  Thermal  Pumping 
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Tubular  Exchanger:  Mechanical  Punaping 
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Figure  15d.  Total  pressure  drop  and  pumping  Figure  15e.  Total  pressure  drop  and  pumping  power  in 

power  in  heat  exchanger  system  and  reactor  heat  exchanger  system  and  reactor  vs  speed  in  pipes 

vs  number  of  exchanger  tubes.  Pressure  and  connecting  reactor  with  exchanger  system, 

pumping  power  gene  rated  by  active  part  of  pipe 
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Tubular  Exchanger:  Mechanical  Pumping 


Tubular  Exchanger:  Mechanical  Pumping 
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Tubular  Exchanger:  Mechanical  Pumping 
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Tubular  Exchanger:  Mechanical  Pumping 
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Tubular  Exchanger:  Mechanical  Punnping 


58 


(9>i)  wniNvyn  jo  ssvw 


X  .  X  ^ 
V  0  «  o 
rq 

«  S?  J}  Q 

i-i  • 

«  o  o 

2  J  .2 


I  ? 

>  J  >  J 

s  s  g« 


(d 

•  H  u 


ui 

s  8  E  J 

§  B  .2  B 
•d  .S  S 
d  5  *<-• 

rt  m  ™  O 

h  «  >5 

O  "  O  « 

“  «  S  2 

«  o>  ” 

«  41  5  n 

i-ii 

<43  OJ  41  (0 
V  jO  00  s 

«J  C  J3 
> 


V  4^ 

d  d 

bo  bc 


INNER  RADIUS  OF  EXCHANGER  TUBES 
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Figure  21  d  (Right).  Mass  of  uranium  in  and  volume  TOTAL  LENGTH  OF  RPES  CONNECTING  REACTOR  WITH  HEAT  EXCHANGER  (CM) 

of  heat  exchanger  system  vs  length  of  pipes  connect¬ 
ing  reactor  with  exchanger.  BNL  Log  No.  D-2081. 
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Tubular  Exchanger:  Mechanical  Pumping 
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Tubular  Exchanger:  Mechanical  Pumping 


Figure  23a.  Mass  of  uranium  in  and  volume  of 
heat  exchanger  proper  vs  inner  radius  of  tubes. 
BNL  Log  No.  D-2085. 
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Figure  23b.  Mass  of  uranium  in  the  volume  of  heat  exchanger 
proper  vs  number  of  tubes  in  exchanger.  BNL  Log  No.  D-2086. 
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umeof  headers  vs  inner  radius  of  tubes.  heat  exchanger.  BNL.  Log  No.  D-2089. 

BNL  Log  No.  D-2087. 


Tubular  Exchanger:  Mechanical  Pumping 


Spray  Exchanger:  Mechanical  Pumping 
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Figure  25a.  Temperature  drop  between  inlet  bismuth 
and  outlet  salt  vs  time  bismuth  drops  fall. 


Spray  Exchanger:  Mechanical  Pumping 
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RADIUS  OF  BISMUTH  DROPS  (CM) 

Figure  26e.  Temperature  drop  between  inlet 

bismuth  and  inlet  salt  in  heat  exchanger  vs  Figure  27a,  Height  of  active  part  of  heat 

time  bismuth  drops  fall.  exchanger  vs  time  bismuth  drops  fall. 


.anger  vs  radius  of  bismuth  drops. 


PARAMETER:  DROP  RADIUS. 


RADIUS  OF  BISMUTH  (CM) 


Spray  Exchanger:  Mechanical  Pumping 


power  in  heat  exchanger  system  and  reactor  vs  time  bismuth  drops  fall, 

time  bismuth  drops  fall. 


Exchanger:  Mechanical  Pumping 
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Spray  Exchanger:  Mechanical  Pumping 


radius  of  bismuth  drops. 


Spray  Exchanger:  Mechanical  Pumping 
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Figure  3 Of.  Total  pressure  drop  and  pumping  power 
in  heat  exchanger  system  and  reactor  vs  speed  in  pipes 
connecting  reactor  with  heat  exchanger. 


spray  Exchanger:  Mechanical  Punnping, 


Figure  31a.  Mass  of  uranium  in  and  volume  of  whole  Figure  31b.  Mass  of  uranium  in  and  volume  of  whole 

heat  exchanger  system  vs  time  bismuth  drops  fall.  heat  exchanger  system  vs  time  bismuth  drops  fall. 

BNL  Log  No.  D-2090.  BNL  Log  No,  D-2091. 


Spray  Exchanger:  Mechanical  Pumping 


re  31c.  Mass  of  uranium  in  and  volume  of  whole 
exchanger  system  vs  time  bismuth  drops  fall. 
BNL  Log  No.  D-209?- 


Spray  Exchanger:  Mechanical  Pumping 
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system  vs  time  bismuth  drops  fall, 

3NL.  Log  No,  D-2095.  Figure  31f.  Mass  of  uranium  in  eind  volume  of  whole 

heat  exchanger  system  vs  time  bismuth  drops  fall. 
BNL  Log  No.  D.2096. 


Spray  Exchanger:  Mechanical  Pumping 
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Spray  Exchanger:  Mechanical  Pumping 
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Spray  Exchanger:  Mechanical  Pumping 


pool  at  bottom  of  heat  exchanger  vs  radius  of  pool  at  bottom  of  heat  exchanger  vs  radius  of 

bismuth  drops.  BNL  Log  No.  D-210Z.  bismuth  drops.  BNL  Log  No.  0-2103. 
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